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^: Abstract 

The neutrino emission due to formation and breaking of Cooper pairs of pro- 



^SJ ' tons in superconducting cores of neutron stars is considered with taking into 

^N . account the electromagnetic couphng of protons to ambient electrons. It is 

Qs ' shown that collective response of electrons to the proton quantum transition 

o' 

^^ ' contributes coherently to the complete interaction with a neutrino field and 

O 
r~| ' enhances the neutrino-pair production. Our calculation shows that the con- 

Q^' tribution of the vector weak current to the uu emissivity of protons is much 

rS^ • larger than that calculated by different authors without taking into account 

[^ • the plasma effects. Partial contribution of the pairing protons to the total 

X 

j^ ■ neutrino radiation from the neutron star core is very sensitive to the critical 

temperatures for the proton and neutron pairing. We show domains of these 

parameters where the neutrino radiation, caused by a singlet-state pairing of 

protons is dominating. 
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I. INTRODUCTION 

When the temperature inside a neutron star core is lower than the critical temperature 
Tc for nucleon pairing, the nucleon matter exhibits a condensate of Cooper pairs, which has 
thermal excitations in the form of not paired quasi-particles. Cooper-pair formation and 
pair-breaking coexist in statistical equilibrium and result in additional neutrino-pair emission 
from the neutron star. Under certain conditions, neutrino emission due to Cooper pairing 
of nucleons may be sufficiently high, compared with or even larger than the emissivity of 
Modified Urea process in non-superfiuid matter. This mechanism of energy loss was proposed 
by Flowers et al. many years ago |l| but have been included into cooling simulations only 
recently |^, P], @]. It was shown, that neutrino emission due to formation and breaking of 
Cooper pairs can accelerate the neutron star cooling. 

Flowers et al. ||l[] have studied the neutrino-pair emissivity for the case of a singlet-state 
pairing of neutrons. One more independent calculation of this process was suggested by 
Voskresensky and Senatorov 0. They have obtained a similar equation for the neutrino- 
pair emissivity, but contrary to previous authors, the expression obtained by Senatorov 
and Voskresensky contains the axial-vector contribution. A new calculation, performed by 
Yakovlev et al. [^], has confirmed that the axial- vector contribution arises actually due to 
relativistic effects and is negligible in the case of singlet-state pairing of neutrons. 

Yakovlev et al. |^ have considered also the more probable case of a triplet-state neutron 
pairing. 

The case of proton pairing was considered in [0. Specifics of the proton pairing, which 
we are going to discuss in this paper, occurs due to a smallness of the weak vector coupling 
of a proton. Cooper pairing of protons take place likely in ^^o-state |§]. When protons 
are treated non-relativistically, the total spin of the Cooper pair in the singlet-state is zero. 
By this reason, the axial-vector contribution of the proton weak current to the neutrino 
emissivity occurs as a relativistic correction, which, however, is found to be larger than the 
vector contribution due to the numerical smallness of the above weak vector constant of 



protons [0]. 

In our paper we argue that, in spite of the smallness of the proton vector weak couphng, 
the latter plays an important role in the neutrino-pair radiation from pairing protons. The 
above inference about negligible contribution of the vector weak coupling was made on the 
basis of calculations which ignored electromagnetic correlations among the charged particles 
in a QED plasma. Actually, protons in the plasma are coupled to ambient electrons via 
the electromagnetic field. The electron vector weak coupling with a neutrino field is much 
stronger than that for protons, therefore the virtual electron cloud, associated to the medium 
polarization, strongly modifies the effective vector weak current of protons 0. The collective 
effects caused by polarization of the electron plasma manifest themselves in the Debye volume 
around the proton. As the wavelength of radiated neutrino and antineutrino is much larger 
than the electron Debye screening distance D^. (Typically, the ratio is of the order of 10 
or larger.), the induced virtual electron excitations generate neutrinos coherently with the 
protons undergoing the quantum transition and leads to stronger neutrino emission than 
that due to the direct proton interaction with a neutrino field. 

In the present article we suggest the calculation of pu radiation caused by Cooper pairing 
of protons in the ^S'o-state, taking into account the above collective effects. We consider the 
medium as consisting of neutrons, protons, and electrons. Generalization including muons 
is also considered. 

To incorporate the electromagnetic correlations, one has to take into account exchange of 
photons between charged particles in the medium. For this we first consider some properties 
of the electromagnetic field in the plasma consisting of superconducting protons and normal 
degenerate electrons. In Sec. |I| we discuss the corresponding electromagnetic field equations, 
polarization functions of the medium, and derive the in-medium photon propagator. Further 
we demonstrate the Lagrangian to be used for weak interaction of protons and electrons 
with a neutrino field. Here we discuss also the problem of nuclear renormalization of the 
weak proton vertex. In Sec. Q we derive the general formula for neutrino-pair emissivity 



by the use of the Optical theorem. In Sec. |3 we apply the closed diagrams to calculate 



the emissivity of the process in the loop approximation reproducing the result obtained by 
different authors, without taking into account the collective effects. In Sec. |V|we calculate the 
neutrino-pair emissivity by the use of the Random phase approximation, which incorporates 
the electromagnetic correlations in the medium. We derive the expression for the RPA 
neutrino-pair emissivity. Some practical formulae and numerical estimates are given in Sec. 
V|. In Sec. [VI I| we show numerical results for some typical parameters of nuclear matter, 



and discuss the efficiency of neutrinos due to proton pairing. Summary and conclusion are 



in Sec. VIII. 



In what follows we use the system of units in which h = c = 1 and the Boltzmann 
constant ks = I- 

II. GENERAL FORMALISM 

To incorporate the collective plasma effects, one should take into account exchange of 
photons between charged particles in the plasma. For this we consider first the electromag- 
netic properties of the medium consisting of superconducting protons and normal degenerate 
electrons. 

A. Electromagnetic field equations 

The superconducting electric current of protons and the normal electric current of elec- 
trons coexist, participating in electromagnetic oscillations of the medium. At zero temper- 
ature, the total Lagrangian density of both the complex field "^ of proton Cooper pairs and 
electrons interacting with the electromagnetic field, has the form 

where k is a constant for proton-proton interaction resulting in the proton pairing, and 
M^ ~ 2M* is the mass of the Cooper pair consisting of two protons of effective mass M*; 
F^u = d^Ai, — dyA^ is the tensor of electromagnetic field; j^ is the electron current, and 



Lg is the Lagrangian density of free electrons. As a result of spontaneous breaking of the 
symmetry of the ground state of the system, the vacuum expectation value of the Cooper 
pair field is nonzero ( |\E'| ) =^o- Since the ground state of the system corresponds to zero 
spin and zero total momentum of the Cooper pair the vacuum expectation value is connected 
with the number density of Cooper pairs Ng by the relation Ns = 2E'$l, where E ^ 2M* 
is the energy of the Cooper pair corresponding to zero total momentum, and A^^ is one-half 
the number density of paired protons Np. Thus, 

With allowance for the interacting with the electromagnetic field A^, we represent the field 
of the Cooper condensate in the form 

^ = ^o(l + p)exp(-2«e0), (3) 

where p and (p are arbitrary real functions of the coordinates and time, and substitute in 
the Lagrangian density ([1|) for the fields p and (p. Taking into account the following identity 

\{d, + 2teA,) ^\' = ^l {d,pf + Ae'^l {A^ - d,<j>f (1 + pf , (4) 

we can make the gauge transformation A'^ = A^ — d^cj). Since the quantity F^u = F'^^ = 
d^Al — duA'^ as well as the electron current j^ = j'^ are gauge-invariant we obtain 

L = ^0 {d,pf + Ae'^>lA^^ (1 + pf + M?^ |V^|^ - A IV^^ - j^i^l - f,^' + ^°- (5) 

As a result, the Goldstone field is absorbed by the gauge transformation. Considering 
this and taking variations of the Lagrangian with respect to the field A'^, in the linear 
approximation we obtain the equation 

a.9^A; + 32vre2v[;2^; = 4vrj-, d^A^ = 0. (6) 

In the absence of the electron current (j^ = Oj, this equation would describe the eigen 
photon modes of mass 



m. 
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327re2^g = JAnNpeyM*. (7) 



This is the well-known Higgs effect. The Goldstone field has been absorbed by the gauge 
transformation. As a consequence of this, the photon has acquired a mass and an additional 
longitudinal polarization. 

However, the polarization of the electron plasma modifies the massive photon spectrum. 
In fact, in the four- momentum representation, the electron current, induced by the electro- 
magnetic field, is 4:njf^ = —nj^jA"^. Here irj^J is the electromagnetic polarization tensor of the 
electron plasma. At finite temperature, one has to take also into account the electric current 
of the normal proton component consisting of non-paired proton quasi-particles. This can be 
described by the proton electromagnetic polarization tensor tvJpJ. One has 47rj^ = —n^^jA"^. 
Thus, the equation for the field (^ takes the form (In what following we omit prime in the 
notation of the field): 

K^A^ - m'^A^ - vr^, {K) A" = 0, K^'A^ = 0, (8) 

where the electromagnetic polarization tensor is: 

vr.. = vri;) + ix^} (9) 



B. Orthogonal basis 

To specify the components of the polarization tensors, we select a basis constructed from 
the following orthogonal four- vectors 

h^ ^ %il, F^ %^, (10) 

where the space-like unit vector n = k/A; is directed along the electromagnetic wave- vector 
k. Thus, the longitudinal basis tensor can be chosen as L^^ = —Fl'^ , with normalization 
LP = 1. The transverse (with respect to k) components of H'''^ have a tensor structure 
proportional to the tensor T^^ = (g^^ — h^h^^ + FF), where g^^ = diag(l, — 1, — 1, — 1) is 
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the signature tensor. This choice of T''^ allows us to describe the two remaining directions 
orthogonal to h and /. Therefore, the transverse basis tensor has normalization T^ = 2. One 
can also check the following orthogonality relations: IpT^^^ = 0, as well as /cjT*^ = 0, and 
KpLP'' = KpTP^" = 0. 
In this basis 

nfP {K) = ni (K) L^^ + tt^ {K) T"'^, (11) 

where the longitudinal polarization function is defined as tti {K) = (1 — uj'^/k'^) tt^^ and the 
transverse polarization function is found to be nt (K) = {gp^nf'^ — ni) /2. 

C. Real part of polarizations 

To the lowest order in the fine structure constant, the polarization tensor of the electron 
gas is defined as follows 

n(-)pp{K)=Ame'TTj^^rG{p)YG{p + K), (12) 

Here G{p) is the in-medium electron Green's function. This is so-called one- loop approxi- 
mation corresponding to a collisionless electron plasma. 

In the case of a strongly degenerate ultrarelativistic electron plasma {vp — 1), one has 



Here the electron plasma frequency and the Debye screening distance are defined as 

2 '^ 2 2 4:7inf,e ^ o 2 f-ic:\ 

"- = ^^ /^e = ^^, ^ = 3..,., (15) 

with fie and Ue being, respectively, the chemical potential and the number density of elec- 
trons. Our Eq. (0) differs from Eq. (A39) of the Ref. 0] by an extra factor [uj'^/k'^ — 1) 



because our basis l'^, h^ is different from that used by Braaten and Segel by the factor 
y/K^/k. All components of the complete tensor (pJ]) identically coincide with that obtained 
in W^ for the ultrarelativistic case. In the case u > k, interesting for us, when the Landau 



damping is forbidden, the one-loop electron polarization functions tc^I are real-valued. 

The real part of polarization functions of the normal proton component is caused by 
dynamics of non-paired proton quasi-particles. For uj > k this contribution is of the order 
ofH: 

„ („) 1 (kVpy 1 kVp , 

where Dp ~ D^, is the Debye screening distance of protons. This proton contribution is 
negligible compared with the contribution of ultrarelativistic electrons. Thus, we neglect 
the real part of proton loops and assume that the real part of the polarization is defined by 
the electron contribution (|l^), (0) only: 

Reni = 7i\'\ Re7r, = 4^) (17) 



D. Imaginary part of polarizations 

The imaginary part of polarizations, responsible for the in-medium photon decay, is 
caused by formation and breaking of the proton Cooper pairs. 




0*0 



Fig 1. Feynman graphs contributing to polarization of the proton plasma. Two proton 
loops, in the right-hand side, are contributions from the normal quasi-particle propagator 



and from "so-called" anomalous propagator which arises in the superconducting phase [12 



The simplest way to calculate the imaginary part of the polarization tensor is an utilizing of 
the unitarity relation for conversion of a photon into itself through intermediate state with 
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two quasi-particles. This corresponds to the loop approximation for the proton polarization 
tensor, as shown in Fig. 1. We consider the retarded polarization tensor. At finite temper- 
ature, the unitarity relation with the intermediate state with two quasi-particles has to be 
written as 

ImTi'"' =Ane^J2^^i^v^*v)^ (18) 



exp (^^ j - 1 ij 

where Aiy is the matrix element of the vector transition current caused by formation of the 
Cooper pair, i. e. due to annihilation of two proton quasi-particles of momentum-spin labels 
(p,s; p', —s) and the total four-momentum K. The summation symbol includes statistical 
averaging and summation over initial and final states of the background, i. e. summation 
over all initial p = {E, p) and final p' = {E', p') states of the proton quasi-particle has to 
take into account the Pauli principle via the appropriate blocking factors, with the Fermi 
distribution function f{E) = l/[exp(i5/T) + 1]. 

The nonrelativistic limit can be obtained by replacing ip'y^ip — > \l/+^, ip'-fiip -^ , 
where ^ is the secondary-quantized nonrelativistic spinor wave function of quasi-protons 
in superconducting matter, which is determined by the Bogoliubov transformation. The 
transformation for singlet-state pairing is of the form (See e.g., ||12|| ): 



^ = ^ exp (zpr) Up (Xsttp,s + X-sOp-s) 
p 

+ J2 exp (-?pr) Vp (xs«p,_, + X-s«p,s) , (19) 

p 

where p is a quasi-particle momentum, and 

E = Ve2 + A2 (20) 

is its energy with respect to the Fermi level. Near the Fermi surface, at \p — Pf\ ^ Pf, one 
has e = Vf{p — Pf)] s enumerates spin states; A is a superfiuid gap at the Fermi surface 
(A ^ pf Vf); Vf = pf/M* is the proton Fermi velocity; Xs is a basic spinor {xtXs' = ^ss')] 
cxp^s and di" ^ are the annihilation and creation operators, respectively. The coherence factors 
are 
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The squared matrix element for the proton-pair vector transition has been calculated by 
Flowers et al. []I| 

J2 I (P,^; p', -a ^+^ O) ' = 2 (Mpv + Up'Vpf . (22) 



Thus, for non-relativistic protons we have 



Y.Tr{M'yM*yn =9"'9"'J^. 



,/ ^ (2vr)3 (2vr)3 

{27i)W^'\p + p'-K) f{E)f{E')2{upVp, + Up,Vp)\ (23) 

The dominant contribution to the integrals comes from proton momenta near the Fermi sur- 
face. As the neutrino-pair momentum is A; ~ Tc ^ Pp, one can neglect k in the momentum 
conservation (5-function, thus obtaining p' = —p. After this simplification, all integrations 
become trivial. Assuming u; > 0, we obtain 

Im.- = ^^ME^tanl. (^) e{u - 2A), (24) 



cuy/(cj2 - 4A2) V4ry 

where 9{x) is the Heaviside step-function, M* is the effective proton mass, and 

Im Ti^ = Im 7r°^' = Im tc'^ = (25) 

By the use of decomposition of the polarization tensor 



n 



PH 



(K) = -ni (K) LP^" + TTt (K) Tf" (26) 



we obtain finally: 






Im7rt(K)=0, (27) 
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E. In-medium photon propagator 

The in-inedium photon propagator, taking into account the medium polarization due 
to superconducting and normal induced currents, can be found by summation of all ring 
polarization diagrams. This can be performed with the aid of the Dyson equation graphically 
depicted in Fig. 2. 








Fig 2. The graphical representation for the Dyson's equation. The in-medium photon, 
shown by thick dashed line, is the infinite sum of ring particle-hole diagrams, connected by 
thin dashed line of the massive photon in the superconductor. 

Here the thin dashed line means the above massive photon, which, in the wave number- 
frequency space has the propagator of the standard form 

D^'i^K) = -——,(g^^---^]. (28) 



m^ \ m^ 



An analytical form of this equation is 



Solution is: 



D^'' (K) = Di^" (K) + ^Di;' (K) TT,, (K) D^'^ (K) . (29) 



4:71 Rf^K^" 
D^'^ (K) = Di (K) L^" + Dt (K) T'^^ - -^—p^, (30) 



where the longitudinal photon propagators is: 



K^ 



^' (^) = 17-2 ^ 1 ■ (31) 

A ^ — m^ — Re TT; — Z Im TT/ 

For uj > k, the imaginary part of the transverse polarization vanishes, therefore the retarded 
transverse propagator has the form: 
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^'W- A-^-,„;-\e.. + .o (^^' 
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Due to conservation and gauge invariance of the vector current IKxtx^^ = tv^^K^ = 
the last term of Eq. ( ^0|) can be omitted. 

According to Eqs. (|3lD , (^2|) , the longitudinal and transverse electromagnetic field obey 
the following equations. 

(-o;2 + e + ni (to, k) + m^) $ (cu, k) = 0, (33) 

(-tu^ + P + vTt (cj, fc) + m^) A {u, k) = 0. (34) 

As RevT; J ~ cUpg ^ m^, the high-frequency polarization of the medium is mostly caused by 
the electron gas, and the superconducting condensate of protons makes only minor correction 
to the eigen photon modes in the medium. The longitudinal eigen mode of oscillations has, 
however, a damping because of absorption of the photon by proton Cooper pairs. 

One can clarify the physical nature of the photon mass by considering a static limit. In 
this case, presence of the superconducting condensate is crucial for the vector electromagnetic 
field. As RevTt (0, k) = 0, Eq. ( ^l)) has the following static form 

(fc2 + m^) A (0, k) = 0. (35) 

This equation has no non-zero solutions for real k. Thus, the static magnetic field in the 
medium must vanish. This is well known the Meissner effect. The field 

produced in the medium by any external static electric current }ext{k), is screened in a 
distance r as exp {—m^r) . Therefore, the photon mass m^ can be identified with the inverse 
penetration depth of the superconductor. 
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F. Weak interactions 

For convenience, we write the low-energy Lagrangian of the proton and electron interac- 
tion with a neutrino field in vacuum as follows 

Here Gp is the Fermi coupling constant, and the neutrino weak current is of the standard 
form 

Jm = ^7m (1 - 75) J^- (38) 

By these notations, the weak current of a bare proton 

j; = ^I'iCv - ^^75)^, (39) 

where ip stands for the proton field, is constructed with reduced constants of the vector 
Cy = 0.5 — 2sin^6'(4/ ~ 0.04 and axial-vector Ca = 0.5(7^ coupling; qa = 1.26, and 9w is 
the Weinberg angle. Reserving the capital letter notations Cy and Ca for proton coupling 
constants, we will use, at the same time, small letters, cy, for electron coupling constants 
with Cy = 0.5 -|- 2 sin^ 9]^ ~ 0.96 for electron neutrinos, and c'y = —0.5 + 2 sin^ 9w — —0.04 
for muon and tau neutrinos. The electron weak current is then of the form 

j^ = uY{cy - CA15)U. (40) 

Strictly speaking, the proton weak coupling vertex should be renormalized due to 
nucleon-nucleon correlations in nuclear matter. However, in vertices with neutral currents 
the nucleon correlations are weakened [^, and are negligible with respect to the plasma 
effects under consideration (See Appendix A). Therefore we will neglect the nuclear renor- 
malization of the weak nucleon vertex. As for the electromagnetic vertex of the proton, its 



renormalization is also negligible for the matter density of interest [^ 
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III. NEUTRINO EMISSIVITY 

We consider the total energy which is emitted into neutrino pairs per unit volume and 
time. By the Optical theorem, for one neutrino flavor, the emissivity is given by the following 
formula: 

where the integration goes over the phase volume of neutrino and antineutrino of the total en- 
ergy u = UJ1 + UJ2 and the total momentum k = ki + k2. In this formula U'^'^ is the retarded 
polarization tensor of the medium which has ends at the weak vertex (Cvlfi ~ C*^7/i75)- 
In the absence of external magnetic fields, the parity-violating mixed axial-vector polar- 
ization does not contribute to the rate of neutrino-pair production. In fact, by inserting 
/ d'^KS'^^^ (K — ki — k2) = 1 in this equation, and making use of the Lenard's integral 

r fhfh^^,) (K-k,- k^) Tr {ff*) 
= ^ {k,K, - K'g,^) e {k') 9 (co) , (42) 

where 6{x) is the Heaviside step function, we can write 

G% r 00 

exp , y 

Since the axial-vector polarization has to be an antisymmetric tensor, its contraction in 
with the symmetric tensor K^Ky — K'^g^j_i, vanishes. Thus, we can take the weak polarization 
tensor as the sum of vector-vector and axial-axial pieces. 

IV. LOOP APPROXIMATION 

To simplest approximation one can write 

^:o, = ClK' + ClTi^L (44) 



Q = ^ [ — r\ — i^ n^'^ (K) (k,k, - K'^g^^ e (k^) e {uj) dud'k. (43) 
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where the polarization tensors Hy^ are given by the same diagrams of Fig. 1, but with 
the ends at weak vertex. Namely, the vector-vector tensor Uy" has ends at the vector weak 
vertex, while the axial-axial tensor U'^ begins and ends at axial vertex. 

Obviously, the imaginary part of the weak polarization tensor can be obtained in the 
same way as the imaginary part of electromagnetic polarizations. The electromagnetic 
polarization differs from the weak vector- vector polarization only by a constant factor 47re^, 
which is included into the definition of the electromagnetic tensor. Thus, 

Imnir = T^Im7r^^ (45) 

4:716'' 

By the use of decomposition of the vector weak polarization tensor 

U^^ (K) = III (K) LP^ + lit (K) TP^ (46) 

we obtain 

Imn. (iT) = -^1^ ^££:^tanh f^V(.. - 2A), 
lmnt(ir)=0, (47) 

As mentioned above, in the case of non-relativistic protons, the total spin of the Cooper 
pair in the ^S'o-state is zero, therefore the axial contribution is proportional to a small param- 
eter Vp. However, due to numerical smallness of the proton weak vector coupling constant 
the axial contribution should be taken into account. Considering the non-relativistic axial- 
vector current of the proton to accuracy Vp we have to include also the relativistic correction 
to its temporal component by replacing ip'^^'^^ip — > \E'+(^ + (^~^\1/, ip'-ji'-^^ip -^ \l/+o"j\l'. Here 
the proton bispinor ip is presented as combination of an upper spinor ^ and lower spinor 
(p = — 'icrjVj\l//(2M), where M is the bare proton mass, and ctj are the Pauli matrices. 

Direct calculation shows that the tensor I'^ = J2a {-^a-^a) is diagonal, with two 
nontrivial elements 0: 

1^2 = If = If = 2{u'v - v'uy. (48) 
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The subsequent calculations with the help of the formula 

2 



lmU^/ = J2T^HM^^M7), (49) 



yield 



exp ( I ) - 1 ij 



MnT^g^Kl f-!^-^^ tanhf^).(.^2A) (50) 



uJ{uj^ - 4A2) V4T 



vr M2 



Imn^i = Imni^ = Imnf = — fcVJ P^ t ^^ tanh (— ) e{uj - 2A) (51) 



37r ^ a;Vcu2 - 4A2 V4T. 

Once the imaginary part of the polarization is calculated, one can evaluate the uu emis- 



sivity as indicated in Eq. (||). By the use of the Eqs. (H), (|4|), (|3), (jSg), (|TD) we 
obtain 

gloop ^ gloop ^ gloop^ (53) 

where the vector weak coupling contribution is 

Q'r = Y^,KGlClpFM*/sJ \ dz, (53) 

ISvrS A (exp (yz) + 1)' ^(^2 _ i) 

with y = A/T and TVjy being the number of neutrino flavors. The axial- vector coupling gives: 

*'-(^^i)''^|«'''°' (^^) 

This result reproduces the emissivity obtained in [|T| and 0] . In spite of the fact that the 
Fermi velocity of protons is small, the axial term ( [5^ ) is typically larger, than the vector 
term (0) due to Cl/C^ > 1. 

V. RANDOM-PHASE APPROXIMATION 

A. RPA polarizations 

In the above loop approximation we ignored electromagnetic correlations among the 
charged particles in a QED plasma. To incorporate the collective plasma effects, one can 
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use the Random phase approximation (RPA). In the case under consideration, the RPA 
means a summation of all ring diagrams connected by the line of a massive photon. One can 
perform the above summation with the aid of the in-medium photon propagator D^^, {K), 
which is solution of the Dyson equation collecting all ring particle-hole diagrams, responsible 
for the electromagnetic polarization of the medium, as shown in Fig. 2. In this way, the 
graphical representation of the RPA weak polarization tensor Yly" is of the form shown in 
Fig. 3. 

<nf^ = 2J^. c5<f> - <3 . C.C.O- - -O 

pee p e 

e P P P 

e p p P 

Fig 3. The graphical representation for the RPA polarization tensor. The in-medium 
photon, shown by thick dashed line, is a solution of the Dyson equation collecting the 
infinite sum of ring particle-hole diagrams, connected by the line of the massive photon in 
superconductor. 

Here the first term describes the contribution, which comes directly from the pairing protons. 
The other terms are caused by the plasma polarization. Analytically one has: 

im nf p^ = ^ [c^ ^^ ^' + i^ (^^ ^^ ^' ^^^ ^')' 

-2CvCv Im ni Re A Re tt; - C^ (Im Tiif Im a)] • (55) 

Here we took into account the fact that 

Re7r;^^ = Im7rf'^ = (56) 

as discussed above. An extra factor (47re^) appears in Eq. (|55D because we replaced weak 
polarization tensors by electromagnetic tensors of the plasma. One more factor (47r) ap- 
pears because the factor An is traditionally included into the definition of electromagnetic 
polarization tensors. The minus sign, in front of terms, which are proportional to CyCv, is 
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due to the fact that those diagrams include the proton charge e > at one of electromag- 
netic vertices of the virtual photon, and the electron charge — e at the opposite end. The 
corresponding factor e^ which appears accordingly to the diagrams with a photon line, is 
also included into the definition of polarization tensors. For the transverse weak polarization 
we obtain 

Im nf PA = ^ Im D, (Re 7^,)^ (57) 



By replacing the photon propagators by their explicit expressions (plD, ( p2|) we obtain: 



imnr^-^"^"' 



H 



Ane"^ 



7^2 - m2 - Re VT/J + (Imyr^)^ 



2CvCv {k^ - m^ - Re vr^) Re vr^ c^ (Re mf 



(^K^-m^-Reni) + (ImTr^)^ (^K'^ - m'^ -Rem) + (Imvr;^ 



(58) 



c2 



Im nf PA = ^ (Re Titf 6 [k^ - m^ - Re n^) (59) 

The imaginary part of the transverse polarization corresponds to the transverse eigen 



photon mode, which has the dispersion defined by Eq. (^J). The eigen-mode frequency 
uJt{k) is larger than the plasma frequency of electrons. Thus, the term Imllj describes 
the contribution from the decay of real transverse photon of energy Ut (k) > ujpe ^ T^. At 
temperatures T < T^. ^^ ujpe-, the number of such photons in the medium is exponentially 
suppressed, therefore this term can be neglected. 

By the same reason we can neglect the RPA corrections to the axial medium polarization. 
Detailed consideration shows that these corrections occur due to exchange of transverse 
photons between charged particles in the plasma. As the imaginary part of the transverse 
photon is the delta-function, the RPA corrections to the axial polarization are of the form 

5 Im nf PA = _£k_ (Re T,^yf S (r^ - m^ - Re vr*) , (60) 

where Re ttav is the mixed axial- vector polarization function of the electron plasma. Obvi- 
ously such a correction is caused by the transverse plasmon decay. Thus, we take 



(lmnT)"" = (lmnr)'°°" (61) 

Further, as the polarization function tci ~ Up^, we can neglect K^ ~ T^ ^ u^^ in Eq. 
]). By replacing also the common factor as specified by Eq. (^) we obtain: 

2 



imnfp^ = imn. 



^2 (^ + Re7r; 



V / X 2 2 

m^+RevT;) + (ImvTi) 



^ K + Re7r,)Re7r, ^ (Revr,)' 



(62) 



[m^ + ReTTij +(Im7ri) [m^ + Reiiij + (ImTr^) 

where the loop polarizations are defined in Eqs. (^T]), p7|), and (|T3|); and m^, given by 
Eq.(|^, is the mass acquired by a photon in the superconductor due to the Higgs effect. 
As noticed above, this mass can be identified with the inverse penetration depth of the 
superconductor. If to eliminate the plasma polarization by taking formally vr^ ^ then 
Eq. (|62D reproduces the result obtained in the loop approximation. However, the plasma 
polarization strongly modifies that result. 

The first term in the brackets of Eq. (162D is the summarized contribution of two diagrams 
of the Fig. 3, which are proportional to Cy. This term is always less than unity. This 
means, that ambient protons screen the weak vector interaction of the pairing protons with 
a neutrino field. 

Let us estimate relative contributions of different terms in Eq. (|62|). According to 
Eqs. (|^ and ([T5|), the penetration depth of the superconductor is typically larger than 
the electron Debye screening distance, i.e. m^ ^ D~^. Thus, m^ < |Re7r/|, and therefore 
m^ + RevT/ ~ |Re7r/|. 

For electron neutrinos one has Cy/cy ^ 1. 74 x 10"^ and Cy/cy'-— 4. 1 7 x 10^^. Thus, for 
electron neutrinos, the first two terms in Eq. ( |62|) are small compared with the term, which 
is proportional to Cy. For muon and tauon neutrinos all the terms are about 10~^ of the 
leading contribution for electron neutrinos. Numerical tests have shown that, to accuracy 
less than one percent, the vector contribution to the emissivity can be described by the 
leading term which is proportional to Cy. Notice that relativistic corrections, which we do 
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not take into account, are of order Vp {Vp ~ 0.06 for a /3-equilibriuni nuclear matter of the 
density p = 2po)- To this accuracy, the above small terms must be ignored. Thus we obtain 
the following analytic expression for the imaginary part of the vector- vector RPA tensor 

Im nr - 4 Im n,- ^^^ -. (63) 

[m'^ + ReTrij + (ImTT,) 

This approximation corresponds to the diagram with two electron loops in Fig. 3. 

B. Including muons 

When the chemical potential of electrons fie is larger than the muon mass m^, the f3- 
equilibrium nuclear matter contains muons. Generalization of Eq. ( |63D to the case of p, e, 
/i plasma is obvious. One has to add the corresponding diagrams with muon loops both to 
the Dyson's equation for the photon propagator and to the RPA polarization tensor. This 
yields 



Y, imfifp^~4imni 



(Revr;^))%2(Revr,.('^)^ 



' J 



2+Re7r;'^)+Re7r;^^) + (ImTrO' 



(64) 



where the muonic polarization function is 

^i'' = ^ f 1 - ^ Vl - ^ In " + f"^ l . (65) 

In Eq. (|6^, the fact is taken into account that the weak vector coupling constant for a 
muon radiating muonic, and tauon neutrinos is the same as that for an electron radiating 
electron neutrinos. By this reason the muonic polarization function comes in the numerator 
of Eq. (|6^) with the factor of 2. The muonic Fermi velocity 



vf = , (66) 



and the Debye screening distance 



H'fj.^Ffi ^ 
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depend on the muonic chemical potential /i^. 

Numerical evaluation under condition of chemical equilibrium yU^ = /Xe shows (see Fig. 
7) that muons introduce a minor correction to neutrino emissivity. 

C. RPA emissivity 

By inserting the imaginary part of the vector- vector RPA polarization tensor 

Im (nf ^) ^' = Im nf P^L^^ (68) 

into the emissivity equation 

Qr = 1^ / .^_^ Im (nr )" {k,k^ - K'g,:) e {k') e (u) dud% (69) 

eXp [rpj i 

we obtain the following vector contribution to emissivity: 

QRPA ^ _GMl r^ rdM^"'"'-'''lmn, <5iljl^ . (70) 

127,^12^ h exp(f)-l (,n? + Re^,) +(Imi,f 

According to above discussion Qy^^., given by this formula, is the uu emissivity of ambient 
electrons perturbed by the proton quantum transition. Therefore it is proportional to the 
square of the electron vector weak coupling constant, and vanishes when the proton energy 
gap is closed, because Imll; (A = 0) = 0. 

Replacing the polarization functions by their explicit expressions (|47|) , (^T]) , and (pTsD we 
finally obtain: 



16 ^2 2 ,.*.7 r , Z-^V^^^l 






[exp(2;y) + 1] 



2 



with 



\x ^f{l-x^)\^{x) ^^^^ 

z^ {z"^ — 1) {I — x^) {ip{x) — \x^) + x^ tanh (zy/2) r^^ 



Wx)^(l-.^)(l-l.ni±|). (72) 
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Additionally to the temperature and the energy gap, the emissivity caused by the vector 
weak coupling depends substantially on two new parameters which does not appear in the 
loop approximation. The parameter 



2 r)2 ^ /ie np Np{T /Tc) 
^ " 3M*n, n, 



A = m'Dl = J^'_^t2PS±l±£l (73) 



"p 



is the squared ratio of the electron Debye screening distance to the penetration depth of the 
superconductor. The latter equality follows from (plSf ) and (0). Here n^ is the number density 
of electrons; Up is the total proton number density, and Np {T/Tc) is the number density of 
paired protons. In the case of p, e-plasma one has np = Ue- Typically A ~ 0.05 -^ 0.1. The 
second parameter is r^ = 2e'^ppM*D'^, which can be written in the following form 

where 

is the inverse Debye screening distance of protons. The parameter rj is typically about one. 
The temperature dependence of the emissivity enters by means of parameter 

A (T) A (0) A (r) 
y = -I-= T. rA(0) (^^^ 



with r = T/Tc. For a singlet-state pairing A (0) /Tc = 1.76 (See e.g. |1T2|), therefore the 
function y depends on the dimensionless temperature r only. Thus, the emissivity (0), in 
[erg s~^ cm~^], can be written as 



or ^ I.IT X lO^H (^) »'. (^) rWfi. 



[exp{zy) + ly 



\2 



where 






M* \noJ ^ ' 
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is the proton Fermi velocity; and uq = 0.16 fm ^ is the number density of saturated nuclear 
matter. 



By adding the contribution (|5^ ) caused by the axial-vector coupling we finally obtain 
the following total emissivity: 

Q = Qf"" + Qa (79) 

with Q^^^ given by Eq. {^ and 

r'y' / , dz (80) 

•^1 [exp(2y) + l]y(z2 - 1) 

VI. NUMERICAL ESTIMATE 

To evaluate emissivities (|77D , (|8y) numerically we use the analytic fit of the energy gap 
suggested in [Q: 

0.157 1.764' 



y (r) = va^7 ( 1.456 - —j^ + —^ , (81) 



The number density of paired protons Np depends on the dimensionless temperature r. 



By the use of the fitted temperature dependence for the energy gap (0) one can calculate 



Np (r) as follows p 



M^) _. , o fdfiE) 



np J dE 



1 + 2 ^^i^de, (82) 



where E = V^+A^, and / (E) = 1/ (exp (E/T) + 1). 
Thus, we obtain: 

/ie Np{T) 



3M* Up ^ ' 



The function Np (r) /np is shown in Fig. 4. 
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Fig 4. Relative fraction of paired protons Np/np versus the dimensionless temperature 
r = T/T, . 

For illustration we consider matter with total baryon density p = 2po consisting on neutrons 
with a small fraction of protons and electrons corresponding to /^-equilibrium. We set the 
effective proton mass M* = 0.7M*. To demonstrate efficiency of the collective effects, the 
Qy^^, given by Eq. ([77|), is plotted in Fig. 5 versus the dimensionless temperature r = T/Tc 
together with that obtained by the loop approximation (LA) without collective effects. One 
can see that, the collective effects substantially enhance, the uu emissivity caused by the 
vector weak current of protons, but the temperature dependence becomes more sharp. Due 
to the collective effects the total emissivity sharply increases when the temperature becomes 
below than the critical temperature for proton pairing. The temperature dependence of 
total emissivity is shown in Fig. 6. 
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Fig 5. Temperature dependence of the RPA vector weak current contribution to neu- 
trino emissivity in comparision with that obtained in the loop approximation (LA). The 
emissivities shown versus the dimensionless temperature r = T/T^ for /?-equihbrium nuclear 
matter. 
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Fig 6. Total emissivity Q = Qy + Qa versus the dimensionless temperature. 



As mentioned above, when the chemical potential of electrons is larger than the muon 
mass, the nuclear matter contains muons. In this case the muonic contribution to the RPA 
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polarization tensor can be taken into account as indicated in Eq. (0). The total emissivities 
calculated with and without the muon contribution are plotted in Fig. 7. The calculation 
was made under the condition for chemical equilibrium /x^ = /ie- 



without muons 




with muons 



0.4 0.5 0.6 0.7 0.8 0.9 1 

T/Tc 

Fig 7. Total neutrino emissivity calculated with and without the muon contribution. 

VII. EFFICIENCY OF NEUTRINO EMISSION DUE TO PROTON PAIRING 

The effect of neutron superfluidity and/or proton superconductivity on different neu- 
trino reactions in the neutron star core is very complicated. Different neutrino production 
mechanisms can dominate at different cooling stages depending on the temperature, and the 
matter density, and vary along with the chosen parameters Tcp, Ten, which are the critical 
temperatures for the proton and neutron pairing. Fig. 8 shows efficiency of different reac- 
tions at the baryon density p = 2po, where the Direct Urea process is forbidden. Neutrino 
emissivities due to the singlet-state proton pairing, and the triplet-state neutron pairing are 
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plotted in logarithmic scale against the temperature together with the total emissivity of two 
branches of the modified Urea processes, and the total bremsstrahlung emissivity caused by 
nn-, np-, and pp-scattering, which are suppressed due to the neutron superfluidity and/or 
proton superconductivity (See [Q, [jT5[, [1^). Two panels of Fig. 8 differ only by different 
choice of parameters T^p and T^n- On the left panel we took T(.p = 5.6 x 10* K, as calculated 
in l]!^ for /?-equilibrium nuclear matter, and Ten = 5.6 x 10^ K, as obtained in. |jl8| for 



triplet-state neutron pairing. On the right panel we assume Tcp = 3.5 x 10^ K - the result 
obtained in |]l9l, and Ten = 8.5 x 10® K, as suggested in pO |. 



O 



ou 
?5 


P = 2po ,-■■'>-'" 




, ''--' /' 1 








/-•'.'' 1 


?0 


/-- .-^ 1 




/■ .^ 






.-/-^ ' 


15 


/7; 


CD CD 1 




o o 




X X 


10 


/ ;■ / 


5 ^ 1 ■ 




LO K? 






II "1 



8.5 9 

Ig T [K] 



9.5 



10 



30 



^ 25 



<f> 20 

O 15 

10 







,' 


P = 2Po 




,'' y' 




,-' 


_ ^ 




x^^^ 






i^ 




^ 


5^ 


■ 


X .. 


^ 




-./ 






/ /. ' 




, , 




^ 


:^ 


' // 


CO 


CT) 


J / ■ 


o 


O 


/ ,' 












/ : 


in 




00 


CO 


/ .' 


II 


Q. 


./ .' 


h-" 


l-y 



8.5 9 

Ig T [K] 



9.5 



10 



Fig 8. Temperature dependence of the neutrino emissivity in different reactions for 
/3-equilibrium nuclear matter of the density p = 2po- The neutrino emissivity due to the 
triplet-state neutron pairing is shown by dashed line. The solid line is the neutrino emissivity 
due to the singlet-state proton pairing. Dot-and-dash line shows the total bremsstrahlung 
emissivity caused by nn-, np-, and pp-scattering; the dotted line exhibits the total emissivity 
of two branches of the modified Urea processes. On the left panel we took Tcp = 5.6 x 10® K, 
and Ten = 5.6 x 10^ K. On the right panel we assume Tcp = 3.5 x 10^ K, and Ten = 8.5 x 10® 
K. 
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Partial contributions of the proton and neutron pairing to the total energy losses are very 
sensitive to the corresponding critical temperatures. Unfortunately the critical temperatures 
Tcp and Ten are not well known up to now. Different authors |]T^ - |Q suggest Tcp and Ten 
which vary in the wide range from 10^ to 10^'' K and sensitively depend on the model of 
strong interactions used for calculation. A large scatter of T^p and T^n does not allow to 
make choice among many different microscopic results for these parameters. Therefore, 
considering the critical temperatures as free parameters, in Fig. 9, we show the regions of 
Tcp and Ten, were different neutrino mechanisms dominate. 
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Fig 9. Regions of Tcp and Ten, in which different neutrino reactions dominate at T = 10^, 
and 6 X 10® K in matter of the density p = 2po 

Additionally to neutrino mechanisms considered above we have included neutrino 
bremsstrahlung due to electron-electron collisions recently calculated by Kaminker and 



Haensel p3[- Unfortunately, their calculation for neutrino emissivity does not take into 
account the collective effects, which substantially modify the effective vector weak current 



of the in- medium electron p^ . However, one can hope that the reappraised result will be 
of the same order of magnitude, and the difference will not be significant in the logarithmic 
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scale we use. 

Two panels of Fig. 9 illustrate the case of the standard cooling at p = 2po for two 
internal temperatures of the neutron star: 10*^ and 6 x 10® K. Our calculations show that 
the represented diagrams does not vary noticeably with the matter density in the region, 
where the proton superconductivity exists. Thus, the diagrams of Fig. 9 show the efficiency 
of the neutrino production due to proton singlet-state pairing in the cooling neutron star. 
One can see that there are wide regions in the diagrams, where the singlet-state proton 
pairing dominates in neutrino production. 

VIII. SUMMARY AND CONCLUSION 

We have calculated the neutrino-pair emissivity caused by a singlet-state Cooper pairing 
of protons by taking into account electromagnetic correlations among the charged particles 
in a QED plasma. To incorporate the electromagnetic correlations we applied the Random 
phase approximation and derived Eq. (^) for contribution of the vector weak current to the 
neutrino emissivity. We found that, virtual electron excitations, electromagnetically induced 
by protons undergoing the quantum transition, generate neutrinos coherently with the initial 
protons and lead to stronger neutrino emission than that calculated by different authors 
without taking into account the plasma effects. Partial contribution of the proton reaction 
to the energy losses from the neutron star is very sensitive to the critical temperatures for 
the proton and neutron pairing. Unfortunately the critical temperatures Tcp and Ten are 
not well known up to now. Their theoretical values sensitively depend on the model of 
strong interactions used for calculation and vary in the wide range from 10® to 10^° K. In 
Fig. 9, for the temperatures interesting for practice, we show the regions of Tcp and Ten, 
where different neutrino mechanisms dominate. The represented diagrams does not vary 
noticeably with the matter density in the region, where the proton superconductivity exists, 
and demonstrate wide domains in the plane Ten, Tcp, where the singlet-state proton pairing 
dominates in the neutrino production. 
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APPENDIX A: NUCLEAR RENORMALIZATION OF THE PROTON WEAK 

VERTEX 

Effect of nucleoli- nucleoli correlations was evaluated in |^. Following this work (Eq. 
2.20), the vector weak form-factor of the proton in nuclear matter is of the following form: 

-PP = (5v - T^f^, (A.l) 

As we use the reduced coupling constants of the proton, the form-factor Upp differs by the 
factor of 1/2 from that obtained in [^]. We have also omited the factors Cq and C^^ which 
finally cancel in the original expression given by the authors. Then, in the actual case of 
k < uj ^ T, the quantity Ann can be written as: 

Ann{uJ,k)=^-(l-'-\n'-±^^ (A.2) 

with s = uj/kVpn- fnp, fnn ~ 1 are constants of the of the theory of finit Fermi systems. 

As the above calculation was made for the case of non-relativistic neutrons we have to 
assume that the Fermi velocity of neutrons Vpn ^ 1- Thus, we have s ^ 1. In this limit 

Anniuj,k)^---<^1 (A.3) 

and the vector weak form-factor of the proton can be simplified as follows 

«pp ^ Cy + ^. (A.4) 



The constants of the theory of finit Fermi systems are not very good known for the (3- 
equilibrium nuclear matter. However, even assuming in the order of magnitude fnp ~ 1 we 
obtain Upp ~ Cy- Thus, the renormalized proton vector weak coupling is found to be much 
smaller than that for electrons. Therefore the nuclear renormalization effects are negligible 
with respect to the plasma effects under consideration. A small nuclear correction to the 
axial- vector weak coupling of protons can also be neglected. 
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